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The trace of integer powers of the local density matrix pv corresponding to the vacuum state 
' reduced to a region V can be formally expressed in terms of a functional integral on a manifold 

with conical singularities. Recently, some progress has been made in explicitly evaluating this type 
P_| ■ of integrals for free fields. However, finding the associated geometric entropy remained in general a 

difficult task involving an analytic continuation in the conical angle. In this paper, we obtain this 
^vjj . analytic continuation explicitly exploiting a relation between the functional integral formulas and 

the Chung-Peschel expressions for pv in terms of correlators. The result is that the entropy is given 
in terms of a functional integral in flat Euclidean space with a cut on V where a specific boundary 
condition is imposed. As an example we get the exact entanglement entropies for massive scalar and 
Dirac free fields in 1 + 1 dimensions in terms of the solutions of a non linear differential equation of 
the Painleve V type. 

1 Introduction 

Recently, there has been a growing interest in different areas of physics for understanding the manifes- 
tations of quantum entanglement in a variety of phenomena. Several measures of entanglement have 
been studied either as a tool to analyze these phenomena or as a subject of investigation in itself. The 
applications range from condensed matter physics [T] and quantum field theory [2] to the physics of black 
holes [3] and holography in the context of string theory [3] . 

In quantum field theory the typical manifestation of entanglement is the statistical correlation of 
£ — , | measurements for independent (commuting) sets of observables located at spatially separated regions. 

A natural measure of the entanglement between a region V with the rest of the space is the geometric 
entropy [5], also known as entanglement entropy. This is defined as the von Neumann entropy S(V) = 
— tr(py logpy) of the reduced density matrix of the vacuum state py — tr_y |0) (0| in the set V (the 
trace in the expression for py is taken over the degrees of freedom lying outside V). 

More generally, there is a one parameter family of information measures also associated to py known 
as Renyi entropies (or alpha-entropies), S n (V) — log(trpy)/(l — n), satisfying the following particular 
limit 

S(V) = \imS n (V) . (1) 

These behave very much like the entanglement entropy. For integer n they are also more suitable to 
explicit calculation, since the traces trpy involved in ([T]), with ne2, can be represented by a functional 
integral on an n-sheeted d dimensional Euclidean space with conical singularities located at the boundary 
of the set V [6] . Here the functional integral is defined on one dimension more than the one corresponding 
to V, which is d — 1. The replicated space is obtained considering n copies of the d-space cut along V, 
and sewing together the upper side of the cut in the k th copy with the lower one of the (k + l) th copy, 
for k — 1, n, and where the copy n + 1 coincides with the first one. The trace of p n is then given by 
the functional integral Z n for the field in this manifold, 

trp" = . (2) 
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In general, solving this integral explicitly is a very difficult problem since we have to deal with a non 
trivial manifold resulting from the replication method. In this scheme, the entanglement entropy follows 
from S n by analytic continuation to n — 1. In the case of free fields, a great simplification follows by 
mapping the n— sheeted problem to an equivalent one in which we deal with n decoupled and multivalued 
free fields [7]. First, we arrange the values of the field in the different copies in a single vector field 
living in a d dimensional space, 

/ 0i (a?) 



(3) 



\ <f>n(x) 



where <f>i(x) is the field on the I copy. Note that in this way the space is simply connected but the 
singularities at the boundaries of V are still there since the vector <3> is not singled valued. In fact, crossing 
V from above or from below, the field gets multiplied by a matrix T or T^ 1 respectively. The matrix T 
is given by 

/ 1 \ 

1 



T 







(4) 



1 

\ (±1)" +1 / 

where the upper sign corresponds to the bosonic case and the lower one to fermions [3 |j. This 
matrix has eigenvalues e l » 27r , with k s = , (n — 1) or kf = — (n — l)/2 , (n — l)/2 depending on 
the bosonic or fermionic character of the field. Then, changing basis by a unitary transformation in the 
replica space, we can diagonalize T, and the problem is reduced to n decoupled fields <pk living on a single 
d dimensional space. These fields are multivalued and live on the euclidean d dimensional space with 
boundary conditions imposed on the d — 1 dimensional set V given by 



f e V . 



(5) 



Here 



and 



are the limits of the field as the variable approaches V from each of its two opposite 



sides in d dimensions. In this formulation we have for fermions 

log(trp^) 



(n-l)/2 

fc=-(n-l)/2 



and for bosons 



log(trp^ 



fc=0 



(G) 



(7) 



where we have written Z[X] the partition function corresponding to a field on a single copy of the 
Euclidean space, which is multiplied by a factor A when the variable crosses V. A related construction 
for two dimensional integrable theories where tip" is expressed in terms of correlators of twisting operators 
is shown in [ID] . 

Explicit calculation of Z [ e l27rfe /™] has been carried out for several cases of interest [7] [HI HI] ■ However, 
the evaluation of the entropy is limited by the difficulties in getting the analytic continuation of the sum 
for non- integer n in order to take the limit ([TJ. 

In this paper we solve this problem by exploiting the relation between ([6|) and ([7|) with the expressions 
obtained by Chung and Peschel [12 for py in terms of correlators for free boson and fermion discrete 
systems. In the continuum limit, these provides the necessary tools to compute the entropy. We find 
that it can also be written directly in terms of a functional integral in flat space. The result is 



S= I dX 



(x-iy 



■logZ[X] 



(8) 



previous work [71 [51 [TT1 we have used the notation Z[a] for the present Z[e* 
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for free fermions, and 

f°° 2 

S = J i dXj^-^logZl-X] , (9) 

for free bosons. Note that the functional integral involved has a boundary condition on the cut V given 
by a real factor (instead of a phase factor as in ^ and (0), being a positive factor for fermions and a 
negative one for bosons, 

(f> !+) (f) = ±X ?&V, A>0. (10) 

In this new scheme, the calculation of the geometric entropy in the set V and in any dimensions is reduced 
to the one of the corresponding functional Z[\], 

In the next Section we derive these equations and, as byproduct, we obtain an interesting relation 
between the d-dimensional free energy log if [A] and the resolvent of a specific correlator acting in the 
d—1 dimensional region V. In Section III we also show how the formulas for the entropy can be obtained 
directly from the expressions (JTJ, © and 0, modulo some assumptions in the behavior of the partition 
function Z[X] as a function of A. An example is given in Section IV where we compute the entanglement 
entropy corresponding to a segment of one dimensional quantum bosonic and fermionic chains which in 
the continuum limit are described by two dimensional scalar and Dirac massive field theories. 



2 Geometric entropy and two point functions for free fields 
2.1 Fermionic fields 

In [12] an expression for py was given in terms of correlators for free boson and fcrmion discrete systems. 
In order to describe the continuum, consider first a set of fermionic fields ^k( x ) which satisfy the canonical 
anticommutation relations 

{*\{x),*M}=8(x-y)& ij . (11) 

The two point correlation functions restricted to V are given by 

(%(x ) y) = (0\*l(x)* j (y)\0) ; x,yeV. (12) 

By definition, the vacuum expectation value (Oy) for any operator Oy localized inside a region V must 
coincide with tr(pyOy), where py is the local density matrix. 

On the other hand, due to the Wick theorem any correlator inside V can be written in terms of 
Cy(x,y). The expectation values calculated with py must also satisfy Wick's theorem. This fixes py to 
be of the form 

py = e- H , (13) 

where H is a local Hamiltonian (modular Hamiltonian) which is quadratic on the fields inside the region 
V, 

H= f dxdyVlWHvfayWy) ; H* t {y,x) = H l3 (x,y) . (14) 
Jv 

The relation tr(pv^\(x) $fj(y)) = Cy(x,y) gives the kernel Hij(x,y) as (we do not write the subscript 
V for notational convenience) 

Hf 3 {x, y) = [log(l - C) - log(C0] y Or, y) ■ (15) 

Then we have 

S - -tr [(1 - C) log(l - C) + Clog(C)] , (16) 

and also 

n — 1 

^~ / pi2Ttk/n \ 

log (trp") = logdet [(1 - C) n + C n ] = £ logdet [C - ei2fffc/n _ 1 j , (17) 
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where this last sum is over the roots of the polynomial (1 — C) n + C n . Recalling eq. (O and comparing 
it with (fTTj) . we conclude that 

logZ[A]=trlog(C-A/(A-l)), (18) 

for any A. 

In order to proceed, we use an integral representation of the logarithm written in terms of the resolvent 

1 1 



logC = - 



d0 



C + + 1 

This gives for the entropy S up to a non universal constant term 

1 1 



d00ti 



C + c-{p + i) 



From equation (|18[) we have 



tr 



1 a dlogZ[A] 
C-A/(A-1) 1 ' dX 



(19) 



(20) 



(21) 



Finally, we use the equation 1(21)1 in (|2"0)l . changing variables — — A/(A — 1), A S (0, 1), in the first term 
of the integrand, and = 1/(A— 1), A S (1, oo), in the second one. These two terms are equal once one 
takes into account that the reflection symmetry of the action on the axis containing V gives place to 

logZ[A] = logZ[l/A] . (22) 

After integrating by parts we get our final expression for the entropy ©. In fact, the corresponding 
boundary term is given by 1/(A — 1) log Z [A] | ^° . It vanishes at A = 1 due to (|22[) (note also that there is 
no contribution to universal terms in S from logZ[l] which is the partition function without a cut). At 
A — > oo it also goes to zero since the free energy logZ[A] increases at most power like in the angle variable 
~ log(A) (we have however no general proof of this behavior). 



2.2 Bosonic fields 

Following similar arguments we have for bosonic fields the following expressions [12 

S = tx 



(C+i)bg(C + i)-(C-i)bg(C-i; 



with C = yXP and X and P the vacuum field and momentum correlators inside V given by 

XV(x,y) = (y)> , 

PV{x,y) = (Mx)n 3 (y)) . 

We also have 



log (trp n ) = -tr 

last sum 
and (|2T))) we conclude that 



1 



n-1 



= ^logdet [C 



1 e 1 



1 



fe=0 



2 P i^l 



(23) 



(24) 
(25) 



(26) 



where the last sum is over the roots of the polynomial (C + |)" — (C 



logZ[A]=trlog(C-i A ' ' 



i) n . Comparing expressions ([7]) 



From 



2A-1, 

and the integral representation for the logarithm we have the entropy in terms of the resolvent 

1 I 



S=- 



d00tr 



C-\+0 C+0+ 



(27) 
lvent 
(28) 



Taking the derivative in (|2T|) we obtain for the resolvent 

f 



tr- 



= {x l)2 dlogZ[X] 



C-(A + 1)/(2(A-1)) v " ^ dX ' (29) 

Finally, we change variables in |28|) as = 1/(1 — A), A G (— oo, 1), in the first term and = A/(l — A), 
A e (0, 1), in the second one, in order to get the entropy ([9|) after an integration by parts. We also have 
to use the reflection symmetry relation 
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3 A different derivation 



We can get the above expressions for the entropy in a different way by the following trick. First we write 
the sum as a contour integral. Consider the following identity 




1 (n-l)f2 / (n-l)/2 \ 

- E w^wtt^ d M £ — T /w= (30) 



(n-l)/2 

1 — n 

fc=-(n— 1)/2 

— i Idzn (V(l/2 - n/2 - nz) - + 1/2 - nz)) /[z] , (31) 

27ri(l — n) J 

where tp is the digamma function and the integral must encircle the interval [—1/2,1/2] on the real 
line. This formula provides an analytic extension of the sum for non-integer n but we have to consider an 
integration contour which includes all the poles of the combination of digamma functions in the integrand 
for non integer n. These are located on the real line and are greater than —1/2. Thus the integration 
contour encircles the [—1/2, oo] segment and we assume that f[z] does not have singularities on it. Besides 
it must decrease fast enough at infinity. Otherwise we can multiply it by a cutoff factor e~ Xz and make 
A — * at the end. We have in the limit n — > 1 + 

i <"^> /2 it./ i i 



Imposing the symmetry f[z] = /[— z] and f[z] ~ z 2 for z around the origin, the pole of the function 
within the brackets for z = do not contribute. Also the function ip'(z) has its poles on the negative real 
axis for z < — 1. Thus we find that only the term proportional to sin _2 (7rz) is relevant for the contour 
integral. Then, if there are no singularities of f[z] for Re[z] > we can choose the integration contour to 
be the imaginary axis. Using ([32|) with f(z) = \og(Z[e l2wz ]) we recover the expression (|SJ) for the fermion 
field entropy. The bosonic case can be treated similarly. 



4 Example: Free massive scalar and Dirac fields in two dimen- 
sions 

In two dimensions, the partition function Z[X] where V is a single interval of length r is determined by 
a non linear differential equation. Following [7] and [5], we introduce functions wd and ws for a Dirac 
fermion and a real scalar, which are related to the respective partition functions Z by 



dr 
d 



r— log Z[e 2 * b }=w D (b,t), (33) 



rj-\ogZ[-e 2 * b ] =w s (b,t), (34) 

where t = mr and m is the field mass. For a fixed b these functions are given in terms of the same 
differential equation for both, the Dirac and the real scalar fields, 



w = -J dy y u (y) , (35) 
1 , u 



-u = 



t l + u 2 

The Dirac and scalar cases differ however on the boundary conditions 



{u r+u{l + u 2 )-j rYT - 2 . (36) 



2 

Ufl(t) — * — smh(bTr) Ki2b(t) as t — > oo (37) 

u D (t) -> -2b(logt+q D ) as t^O (38) 
2 

us(t) — > — cosh(bir)Ki2b(t) as t — > oo (39) 

7T 

u s (t) -> — -— ast^O, (40) 

W t(\ogt + q s ) 1 ^ 



5 



where q D = - log(2) + 2 lE + m±^\ and qs = _ i og ( 2 ) + 2 lE + V-[i/2+ay [i/2-«h] with 7fi the Euler 
constant, ip[a] the digamma function and K the modified Bessel function. The boundary conditions at 
the origin can be deduced from the connection formulas for this Painleve equation [13] or directly from 
the Green function on the cut plane in the massless limit [H [14] . 
The entropic c-function 

c(t)=r±S(r) (41) 

gives the universal part of the entropy S(r) for a single interval in two dimensions [15] . Analogously to 
the Zamolodchikov's C function, it is dimensionless and decreasing with r, and in the conformal limit 
takes the value CV/3, with Cy the Virasoro central charge. We have then 

c D (t) = / db-——w D (b,t), (42) 
J smh(7r6) 2 

cs{t) = / db——-w s (b,t). (43) 
J cosh(7rb) J 

These analytic results can be easily evaluated with high precision by solving the differential equation 
numerically. We have checked they are in perfect accord with previous lattice simulations for c [8j. The 
leading short and long distance terms read 

c D (t) ~ - - it 2 log 2 (i) for t < 1 ; c D (t) ~ ^tKi(2t) for t > 1 , (44) 

CsW ~ I + 2\^(t) for<<<1 ; c s (t) ~ ltK!{2t) fort»l. (45) 

More subleading terms in the short distance expansion for co(i) can be obtained straightforwardly from 
(|42|) and the corresponding series solution of the differential equation [7]. This converges up to t = 1. 

Large distance series expansions in terms of multiple integrals can be obtained from (I42M3| using the 
form factor expansions of the partition function. Remarkably, it has been recently shown that the large 
distance leading term behavior c(r) ~ mrKi(2mr) also holds for a large class of integrable interacting 
theories |10) . 

Other universal terms in different dimensions such as the vertex contribution to the entropy for 
polygonal sets in 2 + 1 [11) , or the leading universal contribution for a long and thin rectangle [5] , can be 
similarly obtained. 
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